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I. INTRODUCTION 



An acoustic ray tracing program is coupled with the Oceanic Boundary' Layer 
Model developed by Garwood [Ref. 1]. The OBL model is a one-dimensional, second- 
order turbulence closure, vertically integrated model of the ocean surface turbulent 
boundary layer, using a two-component turbulent kinetic energy budget with a mean 
turbulent field closure. 

Fisher [Ref 2] investigated the variability and sensitivity of a coupled model 
system. He found that the OBL model, when integrated in time at a single point 
(Ocean Station Papa oO^N 145^\V), predicted inixed-layer structure better than did the 
E.xpanded Ocean Thermal Structure (HOTS) system which was currently in use at the 
Fleet Numerical Oceanography Center (FNOC). McManus [Ref 3] evaluated the 
acoustic performance of a coupled model system at a line of stations in the northeast 
Pacific Ocean. In both cases, the thermodynamic model was initialized with observed 
temperature profiles, and the surface boundary conditions were given by the currently 
available meteorological informations. Then, the thermodynamic forecasts were input 
into acoustic models, such as R,^YMODE or FACT, and the acoustic performance 
was analyzed using the median detection range (MDR) and the convergence zone 
range (CZR). 

This research is the first attempt to link in a single program the OBL model with 
an acoustic model. Simplicity (compared to the operational models available in the 
U.S. Navw) and classification restrictions lead us to develop in Chapter II an algorithm 
for acoustic ray tracing. As no such routine was available at NPS, a copy is attached 
in the appendi.x for further use by students of the Air-Ocean Sciences Department. This 
simple subroutine allows the influence of the atmospheric forcing on the underwater 
sound propagation to be qualitatively analyzed. A summary of the leading principles 
and equations of the OBL model is given in Chapter III. Chapter IV develops in detail 
the actual coupling of the two models into a single computer code. Chapter V gives an 
example of the use of this coupled model applied to a specific area in the western 
.Mediterranean Sea, for diflerent periods of the year having particular acoustic 
properties. .-Ml the simulations were integrated in time using climatological data over 
short time periods varying from a few hours to three days. 
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II. ACOUSTIC MODEL 



A. INTRODUCTION 

Since we assumed in this research that the ocean is horizontally stratified, the 
temperature T and the salinity S are only functions of the depth y and cannot vary 
with the range z. Hence, the speed of sound c is only a function of depth y. In that 
case, according to Ziomek [Ref 4; p.236], the general form of the equation for the 
horizontal range travelled by an acoustic ray is the following : 



(Vq.Zq) are the coordinates of the source, and P(yg) is the initial angle of propagation. 

Thus, theoretically, knowing the sound-speed profile versus the depth, we can 
plot the curve giving the path of an acoustic ray. 

We chose not to use equation (2.1) for the following reasons: 

.According to Snell's law, 

sinpry) _ sinP(yp )_ 
c(y) c(yg) 

and. at a turning point. P(>\p) = tt'2 and c(yjp) = 1/b such that the denominator of 
the integrand of (2.1) goes to zero and the integration cannot be carried out. Thus 
(2.1) is only valid between two turning points. 

Numerically, the integration of (2.1) would be carried out between all the turning 
points of the ray by discretizing the sound-speed profile c(y) with as small a depth 
increment as required to get an acceptable result. Furthermore, during the integration, 
we would have to test for the occurence of a turning point. 




( 2 . 1 ) 



where b is the ray parameter and is given by 



b = sinp(yg) / c(yg). 
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However, a sound-speed profile can be approximated by straight line segments 
matching the profile as best as possible. The more segments chosen, the better the 
computation. Roughly, we often choose 10 to 15 segments depending on the shape of 
the profile. Then, as explained in the next section, the integration of (2.1) can be 
solved analytically for each segment of straight line and leads to very simple equations. 




We shall keep in mind that this acoustic model has to be coupled with an 
Oceanic Planetary Boundary Layer (OPBL) model [Ref. 1] which resolve.s the 
temperature, salinity, and depth of the upper mi.xed layer of the ocean using a one- 
meter-step increment for the depth. 

These are the reasons why we chose to discretize the sound-speed profile with a 
one-meter-depth increment, assume the profile to be a straight line segment within each 
depth increment, and we use the equations derived in the next section instead of 
numerically evaluating the integral of (2.1). We found this method easier for handling 
the turning point problem. 

B. SPEED OF SOUND AS A FUNCTION OF DEPTH WITH CONSTANT 

GRADIENT 

The sound-speed profile is given by 



c(y) = Cq + g X (y-y^) 
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(2.4) 



u 

where g is a constant (with units of s > referred to as the gradient since 



dc(y) 

Starting with Snell's law (2.3), differentiation leads to; 

cosp(y) dp = b dc(y) = bg dy. 

Referring to Figure 2.2, it can be easily seen that : 
dv 

-f- = cosp(y). 



(2.5) 



(2.6) 



(2.7) 




Figure 2.2 An infinitesimal element of arc length ds 
at an arbitrary point P along a ray path in YZ plane. 

Substituting (2.7) in (2.6) yields 



tip(y ) 

ds 



= bg. 



( 2 . 8 ) 
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Equation (2.S) indicates that the curvature along a ray path is constant. Thus, the ray 
path is an arc of a circle. From Figure 2.2, we have 



dz = tanp)(y) dy = sin{)(y) 



dy 

cosp(y) 



(2.9) 



Using (2.6), we get 

dz = sinP(v) . (2.10) 

bg 

Integration leads to 



z 



,.p sinP 

•’po ^ 



dp 



( 2 . 11 ) 



where Pg = P(yo) initial angle of propagation at the source, and P = p(y) is the 

angle of propagation at position P. This leads to 



^ ^0 + (cospg - cosp) . (2.12) 

Solving for the depth y in (2.4) yields 



y = To + “7 (c(y) - Cq) , (2.13) 

and. using Snell's law (2.3), we get 

We can verify that (2.12) and (2.14) are the parametric equations (parameter P) of the 
circle given by : 
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(2.15) 




centered at { >q - sinpQ bg , -{Zq 4- cosPq bg) ) and having a radius of 1 bg. 

C. ALGORITHM 

For this research, we used a modified version of the basic FORI R/W subroutine 
R.AY given in appendix A. In this section we arc going to analyze in detail the 
dilfercnt parts of the algorithm used in the subroutine R/\Y. 

The goal of this subroutine is to plot the acoustic ray path, that is, the curve 
giving the range z in kilometers versus the depth y in meters of an acoustic ray 
emanating from a source at a depth \'q at a given initial angle of propagation. '1 he ray 
path plot is presented beside the graph of the sound-speed profile c(y) in m sec versus 
depth y in meters. 

As we mentioned previously, we used a one-metcr-step depth increment and kept 
track of the depth all along the ray by using an increasing or decreasing index k 
depending on the ray going upward or downward. The index k is, in fact, the integer 
value of the depth. 

The angles of propagation P and Pg are referenced from, the Y axis. For plotting 
purposes, we defined the depth vector { yj, i = O.NN } and the range vector { Zj, i = 
0,NX }. At the i’’^ point of the ray path, the depth index k is 



the "initial" angle of propagation is pg, the sound-speed is Cj^ and. its gradient, gj,. is 
given by 



propagation is p, the sound-speed is Cj^^ j and the gradient is gj^ | . An illustration 
of these difTerent parameters is given in Figure 2.3. 




(2.16) 



Sk " ^k+ 1 • ^k • 



(2.17) 



At the (i-i-1)^^ point on the ray path, the depth index is k±l the "final" angle of 
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For that case, according to (2.3) and (2.12), we would have 



>'i+ 1 



= k - 1 



and 



(2.18) 



^ (cospQ-cosP) 

K“ 1 

with 



(2.19) 



‘'k-1 ^source 



( 2 . 20 ) 




Figure 2.3 Parameters describing a one-meter-dcpth increment 
positive gradient sound-speed profile 
for an upward travelling acoustic ray path. 
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Now we are going to examine in detail the dillerent cases occuring in the 
subroutine leading to some IF. ..THEN. ..ELSE... statements, depending on the sign of 
the gradient (or the curvature) and the upward or downward direction of the ray at the 
i^^ point. 

1. Case I : g|., j > 0, gj. > 0. 

At the i^^ point of the ray path Pj, corresponding to the depth index k, the 
gradient above is (see Figure 2.4) 

8k-I ^k‘ ‘^k-l (2.21) 

and the gradient below is 

8k = ^k+1 -"^k- (2.22) 




Figure 2.4 Case 1 : gj^.j >0, gj^> 0. 

We have three different cases to consider depending on the value of Pq with 
regard to m2 and P^. The critical angle P^ is the value of the initial angle of 
propagation at Pj leading to a final angle P = k'2 one meter deeper, that is, leading to 
a turning point one meter deeper. 
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Accordins to Snell's law. 



sinP^ = 



^^ 1 ;+ 1 



(2.23) 



or 



P„ = arcsin . {2.24) 

Let us examine the three cases Pq> tt, 2, Pg < P^ and P^- Pq since a strict 
equality does not apply for real numbers in FORTRAN language. 
a. Pg > 3T, 2 : 

According to (2.3) and since the ARCSIN function gives a result between 0 
and 71 '2, the angle of propagation at j is : 

P = 7t - arcsin(bcj._|) . (2.25) 

Substituing for P in (2.12), we get the coordinates of the next point Pj 4 - j : 

yi+l = k-l ( 2 . 26 ) 

and 



Z-4.J = Z; + (cosPg + cos[arcsin(bcj^. j)]} 



with 



b = sinp,^,,,.^^ , c 



source ' '"source 



(2.27) 



(2.2S) 
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have : 



pQ Pc • 

This time, using the same equations (2.3) and (2.12) at the point ^^ 4 - 1 . 



P = arcsin(bc^^ j) , 



( 2 . 29 ) 



>’i+ 1 



= k+1 



(2.30) 



and 



4 + 1 



= 7 - 



bgk 



{cosPg - cos[arcsin(bcj^^ j)]} 



(2.31) 



Pc PQ ■ 

in this case, we have to deal with a turning point within the segment 

Pi^i+l' 

Let us derive the general formula for computing j when a turning point 
is encountered. This formula will be used when we study the cases generated by the 
other sign possibilities of the gradient g. 

Applying (2.12) first from Pq to P'. then from P' to P we have (see Figure 

2.5): 



Z' = 2„ H (cOsPr, - COST! 2) 

and 

z = z' + — ^ — [cost: 2 - cos(;t - P„)] , 
b2 



(2.32) 



(2.33) 



or 
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(2.34) 



and 

1 n 

Z = Z -) ; — COS Pn . 

bg ^ 

Substituing (2.34) in (2.35) yields 



(2.35) 



(2.36) 




Figure 2.5 Turning point treatment. 

Going back to our case of interest and using (2.36), we get the coordinates 



■i-H 1 



= y: = k 



(2.37) 
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and 



+ 



bgk 



cospQ . 



(2.38) 



2. Case 2 : g^-i < ® 

This case is ver>’ sinailar to case 1, but now the critical angle is given by : 



p^ = Jt - arcsin(Ck/Ck.i) • 



(2.39) 



Following the notation of Figure 2.6 and applying the same basic equations (2.3), 
(2.12), and (2.36), we have three new cases to examine. 




Figure 2.6 Case 2 ; 

a. Pq < n/2 : 



P = arcsin(bck+ |) , 



(2.40) 
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(2.42) 



z; a. ] = z; • r — 
* * ' bci. 



(cospQ - cos[arcsin(bcj, ^ j)]} . 



b. Pn > Pr : 



P = rc - arcsin(bcj^_j) , 



(2.43) 



1 - 



= k-1 : 



(2.44) 



and 



Zjo.} = Zj + — ^ {cospg + cos[arcsin(bcj..j)]} 
=k- 1 



c. 7T 2 < pQ < P^ ; 



t = = k 

" 1^1 



and 



(2.45) 



(2.46) 



'i+1 



— Z; + 



bgi- 



cosp. 



k-1 



(2.47) 



3. Case 3 : gj^. j > 0» gj^ < 0 

According to Figure 2.7, wc now have only two separate cases to consider, 
depending on the value of Pg with regard to rt, 2. Returning to the two previous cases, 
we will find some similarities. 



a. > Tt/2 : 

This case is similar to case l.a, and equations (2.25), (2.26), and (2.27) 

apply. 

b. Po < jt/2 : 

This case is similar to case 2. a, and equations (2.40), (2.41), and (2.42) 

apply. 




Figure 2.7 Case 3 : gj^_| > 0, gj^ < 0. 



4. Case 4 : <0, 

According to Figure 2.8, we now have four dilTerent cases to examine with 
two possibilities for P^ . 

P = Ti - arcsin (2.48) 

‘^k-l 



or 



P^' = arcsin 



^k+ 1 



(2.49) 
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The four different cases shown above reduce to cases previously analyzed in 
sections for case 1 and case 2. 




Figure 2.8 Case 4 : gj^_j <0, gj^> 0. 



^ Pc' 

This case is similar to case 2.b, and equations (2.43), (2.44), and (2.45) 

apply. 

b. Pc ^ Po "" ^^'2 : 

This case is similar to case 2.c, and equations (2.46) and (2.47) apply. 

Po Pc' • 

This case is similar to case l.b, and equations (2.29), (2.30), and (2.31) 

apply. 

d. jt/2 > Pq > p^. : 

This case corresponds to case l.c, and equations (2.37) and (2.38) apply. 
This last case ends the cascade of IF. ..THEN. ..ELSE... statements that we used to 
write this subroutine. 
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5. Surface and bottom reflections 

The last point we need to discuss is how to handle surface and bottom 
reflections, that is, when the depth index k reaches the values of 0 or N. The value N 
corresponds to the maximum depth given in the profile and is assumed to be the depth 
of the sea floor. In the case where the sound-speed profile does not extend to the 
bottom, it is always possible to extrapolate the sound-speed profile to the bottom using 
the gradient of the profile corresponding to the last straight line segment or, using an 
average value given by the climatology. 





Figure 2.9 Perfect surface reflection possibilities. 
a. Surface reflection 

In this simple acoustic model, we assume a perfect surface reflection. Thus, 
when the depth index k reaches the value 0, we just have to maintain symmetry (i.e., 
the angle of reflection is equal to the angle of incidence) with regard to the horizontal 
to get the correct angle Pg before applying the same equations as in cases l.b, l.c and 
2. a of Figure 2.9 . 
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b. Bottom rejlection 

Following the same idea, we assume a perfect bottom reflection. Thus 
when the depth index k reaches its maximum value N, we maintain symmetry (i.e., the 
angle of reflection is equal to the angle of incidence) with regard to the horizontal (see 
Figure 2.10). 

Notice that the fictitious gradients g.| and g>j have been created in the 
program to handle the cascade of IF.. .THEN. ..ELSE... statements. 




See cases 2.1 and I."!) 




Figure 2.10 Perfect bottom reflection possibilities. 



D. SUMMARY AND DIRECTIVES TO USE THE SUBROUTINE RAY 
1. Main characteristics 

'fhe computation in the subroutine RAY is based on representing an arbitrary 
sound-speed profile by straight line segments in one-meter depth increments and thus 
allows the use of any complicated sound-speed profile to obtain a precise and smooth 
acoustic ray trace. 
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A plot of the sound-speed profile is provided at the same time as the plot of 
the bundle of rays. 

The main assumptions underlying this program are the following : 

• perfect surface reflection, 

o perfect bottom rellection, 

• Hat bottom. 

2. Features of this subroutine 

This subroutine uses many DISSPLA graphics statements. Mence, the main 
program calling subroutine l^AY has to be e.Kecuted using the command DISSPLA. y\ 
choice of multiple display devices is provided to the user through the use of comment 
cards in front of the statements CALL CO.MPRS, CALL TEK61S, or CALL 
CX41(4107). 

The variables required when using R.AY are the following : 

• YO : depth of the source in meters. 

• M : number of ray path desired, 

• BET : array (size M) of the initial angles of propagation Pjource degrees. 

• BO : initial angle of the upper ray of the bundle, 

® DB : increment of initial angles in BET, 

• MM + 1 : number of points in the provided sound-speed profile, 

• CC : provided values of sound speed in m sec, 

• YY : corresponding depth in meters. 

• FLANGE : maximum range desired in kilometers, 

• XN : index of range, 

• Y : array (size NX) of depth in meters, 

• Z : array (size XX) of horizontal distances in kilometers, 

• X ; integer value of the depth of the sea floor in meters. 

• C : array (size X) of sound speed (m sec), 

• G ; array (size X) of its gradient (sec'^), 

• YC : array (size X) of depths (meters). 

The bundle of rays to be traced is defined by the number of rays M, the initial 
angle BO of the upper ray of the bundle , and the angle increment DB between two 
rays. Inside the subroutine, we can also set different equations to define a bundle of 
rays. And finally, by deleting these few lines involved with these computations, we can 
provide our own array BET of initial angles Pspyj-^e- 
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No relation has been derived between NN and the maximum range R.ANGE 
since this relation depends on many dilTerenl parameters such as the Pjource 
sound-speed profile, and the maximum range. If during a plot a ray ends before 
reaching the right side of the graph (especially the steepest rays), a higher number NN 
has to be set. For instance, a value N'N = 5000 seems to handle most of the reasonable 
steep P<;ource ^ maximum range of 30 km. 

The arrays Y and Z have been created for computing and plotting the ray 
paths and the arrays C, G and YC are for plotting the sound-speed profile. 

3. E.xample 

A result of the RAY subroutine is displayed in Figure 2.11. The provided 
parameters and sound-speed profile are listed in Table 1 . 



TABLE 1 

PAR.'\.METERS AND DATA USED IN THE EXA.MPLE 
OF SUBROUTINE RAY SHOWN IN FIGURE 2. 1 1 



^'0 


30. 


YY 0. CC 


1482. 


.M 


10 


30. 


1482.9 


NN 


5000 


55. 


1484.1 


IG\NGE 


30. 


90. 


1480.5 


N 


300 


150. 


1482.0 


MM 


6 


225. 


1484.1 






300. 


1486.2 
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Figure 2.11 Ray tracing provided by subroutine Rv\Y. 
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III. OCEAN MIXED LAYER MODEL 



A. INTRODUCTION 

The study of the oceanic turbulent boundary' layer is a relatively recent field in 
Physical Oceanography. The model used for this research has been developed by 
Garwood [Refs. 1,5). These papers give the formulation of a unified mathematical 
model of the one-dimensional, non-stationary, oceanic turbulent boundary layer. 

The study of these top few tens of meters of the ocean is of considerable 
scientific interest. It influences and can be related to the general circulation of the 
ocean [Ref 6]. The thermal structure associated with this boundary layer must be 
considered when making medium and long-range weather forecast, since a large part of 
the atmospheric energy supply comes from the heat exchanged with the ocean. This 
layer is also a region of primary' biological productivity, which is of significant 
ecological and economic importance. Finally, as an important military application, this 
study can be used in the modeling of acoustic propagation in the ocean, which is the 
goal of this research. 

1. Characteristics of the oceanic mixed layer 

The oceanic mixed layer is that fully turbulent region of the upper ocean that 
IS bounded above by the air-sea interface and below by a dynamicaly stable water 
mass. The wind and intermittent upward surface buoyancy flux through the surface 
(surface cooling at night or in winter time for example) are the sources of mechanical 
energy for the generation of these turbulences. Figure 3.1 gives a general picture 
idealizing density and mean velocity profiles of the ocean iriixed layer. ^ 

2 . Generalities on the dynamics of the mixing 

The depth of the ocean wind-mixed surface layer is typically on the order of 
ten to one hundred meters. The horizontal scale size is that of the internal Rossby 
radius, typically 20 to 50 km. These two dominant scale sizes are much smaller than 
the horizontal scale size of the driving meteorological disturbances, water mass 
features, and distances to lateral boundaries. The approximation of local horizontal 

’in this chapter, in order to be consistent with the notation used in [Refs. 1.5] as 
in most of the geophysical sciences publications, we chose z to be the upward vertical 
axis and (x,y) to be the horizontal coordinates as depicted in Figure 3.1. The vertical 
component of the fluid velocity will be w and the horizontal components will be u and 

V. 
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homogeneity for all mean variables is usually accurate and is a basic assumption in this 
model. 




Figure 3.1 Idealized model for ocean ntixed layer.. 

A sharp density discontinuity of thickness 5 (see Figure 3.1) separates the 
layer from a stable non turbulent thermocline. Minimal stress at the bottom together 
with high turbulence intensity leads also to an appro.ximate vertical uniformity in mean 
velocity and density. We shall note that only small gradients in these mean variables 
give rise to large turbulent fluxes. 

The mechanical energy budget for the ocean mixed layer is depicted in Figure 
3.2. Deepening of the mixed layer is accomplished by entrainment of the more dense 
underlying water into the turbulent region above. This process leads to a potential 
energy increase and cannot take place without an energy source; the turbulent kinetic 
energy of the mixed layer above. 

Retreat occurs when the vertical component of the turbulence is insufficient to 
transport heat, momentum, and turbulence to an earlier-established depth of mixing. 
This happens when downward heat flux (buoyant damping) and dissipation effects 
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Figure 3.2 Mechanical energy budget for the ocean mixed layer. 
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exceed the wind-stirring shear-production of turbulence, creating a stratification of the 
upper ocean. 

Thermal energy and mechanical energy received from the atmosphere not only 
control the local dynamics, but the layer itself modulates the flux of this energy to the 
deeper water masses. Entrainment also converts some of the mean flow energy into 
turbulent energy, over and above the wind-stress production. 

Finally, substantial barotropic and baroclinic features, such as tidal motion 
and internal waves, can be linearly superimposed. The mean fields of concern are 
therefore the horizontally homogeneous components of the total fields. 

3. The model and its features not previously demonstrated 

The vertical and horizontal components of the turbulent kinetic energy (TKE) 
are deternrined implicitely, along with layer depth, mean momentum, and mean 
buoyancy. Layer growth and retreat are predicted. 

Specific features of the TKE budget include mean turbulent field modeling of 
the dissipation term, the energy redistribution term, and the term for the convergence 
of buoyancy flux at the stable interface as shown in a following section. Then an 
entrainment hypothesis dependent upon the relative distribution of the TKE between 
horizontal and vertical components permits the closure of the system of equations. 

The model differs from earlier models in the following ways. First, the 
amount of wind generated TKE to be used in mixing is a function of the ratio of the 
mixed layer depth to the Obukhov mi.xing length L. Second, viscous dissipation is 
dependent on a local Rossby number. Finally, separate vertical and horizontal 
equations for TKE are used, permitting a more consistent interpretation of both 
entraining and retreating mixed layers. 

B. FORMULATION OF THE EQUATIONS USED IN THIS OBL MODEL 
1. Generalities 

The underlying principles employed in studying the inixed layer are the 
combined conservation of mass, momentum, thermal energy, and mechanical energy. 
Most of the physics behind such a one-dimensional model is based on the flux form of 
the Navier-Stokes equations with the Boussinesq approximation. The horizontal 
homogeneity mentioned previously permits the neglect of the horizontal gradients of 
the mean fields. 

Conservation of buoyancy is employed as a generalization of the conservation 
of heat alone. The buoyancy equation is generated from the heat and salt equations 
together with an equation of state. 



p = pf, [1 - a{O-0g) + P(s-Sq)] . 



(3.1) 



and the definition for buoyancy, 



b = giPo-P) Po - 



(3.2) 



where Pq is a representative density at the time and location of consideration, g is 
acceleration due to gravity, a and p are the expansion coefficients for heat and salt. 
All variables are separated into mean and fluctuating components : 

• temperature 0 = T + T' 

« salinity s = S + s' 

» pressure p = P + p' 

• velocities u = L'„ + U + u' 

V = V® + V + V' 
w = \^ 4- w' 



• buoyancy b = B + b' 

Subscript g denotes geostrophic components. 

2. Mean buoyancy and momentum equations 

The first law of thermodynamics for an incompressible fluid and the 
conservation of salt mass, neglecting molecular fluxes, lead to the mean buoyancy 
equation ; 



d^,di = -d'o'Wdz 4- agQ / PgCp . (3.3) 

By invoking the Boussinesq approximation, dropping the negligible viscous terms, 
assuming incompressibility, and subtracting the geostrophic equations from the total 
momentum equations, we obtain the mean momentum equations ; 



d\: dt = IV - du’v.\ dz (3.4) 



8\':ci = -ITJ - dv'w'i'dz . 



(3.5) 
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Integration of (3.3), (3.4) and (3.5) over the entrainment zone from z=-h-5 to z = -h 
leads to the jump conditions for the turbulent fluxes at the bottom of the mixed layer : 



-b'w'(-h) = AB 5h at (3.6) 



-u'w'(-h) = AL' ail at (3.7) 



-v'w'(-h) = AV ah, at (3.S) 

where 

b’ = agT' - Pgs' (3.9) 

and 

AB = agAT - pgAS . (3.10) 

The assumption of vertical homogeneity in the mixed layer permits the integration of 
(3.3), (3.4) and (3.5) from z = -h-6 to z = 0, including the effects of entrainment stresses 
(3.7) and (3.S), and entrainment buoyancy flux (3.6). This yields the bulk relationships 
for mean buoyancy : 

h a<B> at + AAB ah, at = agQ^ / pQCp - ^( 0 ) (3.11) 



ha<u>,at + AAu ah at = -fh<v> - u'w'(O) (3.12) 



ha<v> at + AAV ah at = -ni<u> -v'w'(O) (3.13) 
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where 






(3.14) 



The value of Qg is the net solar irradiance at the surface, minus the long wave back 
radiation, minus the sensible heat flux, minus the latent heat flux. 

The function A is the Heaviside step function : 

A = 1 if ch ct > 0 , 

A = 0 if dh.ci < 0 . 

The brackets < > denotes a vertical mean through the mixed layer, 



As the time step used in the m.odel is one hour, the surface boundary conditions are 
prescribed hourly ; 



< > 



h+O -u-u 




(3.15) 



and denotes horizontal mean. 




(3.16) 



-U'W'(O) = T,^(t) pg 



(3.17) 



-V'W'(O) = Ty(t),Pg 



(3.18) 



-b'w'(O) = g ( ps'w'(0,t) - awT'(0,t) ] 



(3.19) 



with 



^ = P.-,c,u 



a^d^lO 



(3.20) 
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where p,^ is the air density, a drag coelTicient and Ujq the wind speed at 10 meters, 
and with 



s'w'(O) = S(E - P) 



(3.21) 



where E-P = evaporation minus precipitation in cm sec. 

3. Mean turbulent kinetic energy equation 

Subtracting the scalar product of (u,v.w) with the mean momentum equations 
from the mean equation for the total mechanical energy yields the mean TKE 
equation; 



IcE 
2 dt 



- U W :: h 

oz 







0 



where (I) 



(ID (III) (IV) 



(3.22) 



E = 



u'^ - 



V 



W " . 



(3.23) 



The budgets for the individual components of TKE can also be formed : 



1 ^u'' 

2 Cl 



—r~,dU d \v'u' 



p'cu' c 
‘ 3 



n,u'v' - f2.,u'w' 
2 



1 dv'~ — — cV d \v'v'~ , pev' 

_ _ = .V w ~ T“ 

2 ct oz cz 2 PqCv 





2 Cl 



. w 



w p^ 



b'w' - — (— + — ) 

cz 2 Po 



p’cw' c 
PqCz ■ 3 



+ 



ri,u'w' 



(3.24) 



(3.25) 

(3.26) 



(V) (IV) 

where 12 = (12^.12.,, 12^) is the rotation vector for the earth. 

The time rate of change of TKE is usually smaller than the other terms and 
may be neglected. The term (I) represents the rate of mechanical production and is the 
donrinant source of TKE. It is the rate of conversion of mean to TKE by the down- 
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gradient turbulent flux of momentum. The term (II) represents the buoyancy Ilux and 
can be either a source or a sink. This term can become an important source as in the 
case of strong convective cooling in the autunm. The term (III) is the divergence of 
the turbulent Ilux of TKE or turbulent diffusion of itself Locally, at the bottom of the 
layer during occasions of entrainment, a net convergence of flux of energy is necessarx’ 
to maintain the downward buoyancy flux for deepening the mixed layer. The term (IV) 
represents viscous dissipation. Because local isotropy is assumed for the dissipation 
range, £ is divided equally among the component budgets. The terms (V) are very’ 
important terms which sum to zero by continuity in the TKE equation. They cause a 
redistribution of energy among u'', v'* and w'-. The terms (VI) are also redistribution 
terms, but they are due to rotation of the earth. In this study, we will neglect them 
because of the usually short integral time scale in comparison with the time periods of 
interest. 

Assumption of vertical homogeneity in the mixed layer permits the vertical 
integration from z = -h-6 to z = 0 of (3.24) + (3.25) and (3.26). By this step we 
introduce a new variable h. the depth of the mixed layer. 

Up to that point, if we set q". the horizontal component of TKE, such as 

q‘ = u ■ + V *• . (3.2/) 

wc have six variables h, <B>, <U>, <V>, <w'"> and <q^> for five equations 
(3.11), (3.12), (3.13), vertical integral of (3.24) + (3.25) and (3.26). Therefore a sixth 
equation is needed to close the system. Besides, we also need a suitable modeling of 
the different terms of the integrated TKE equation. Following Garwood's arguments 
[Refs. 1,5]. a synthesis of these derivations is given in the next sections. 

4. Modeling of the different terms of the integrated TKE equation 
a. Shear production and turbulent diffusion 

The vertical integral of terms (I) and (III) may be combined to give the net 
"wind-generated" rate of production : 

G = J_^_g[(I) + (IlI)]dz (3.2S) 



or 
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G = m 3 U'’-^ + 1 2 [{AU)‘ - (l\V)-] oh ot (3.29) 

where ir is the friction velocity defined as : 

u* = [u'v'(0)“ -T- v'w'(0)“]''^ . (3.30) 



b. Net buoyant damping 

The integral of (II) over the mixed layer gives the net buoyant damping for 
the whole layer : 




(3.31) 



or 



B = ■ b'W(O)] 



fO .Qh 

piic; J-h-s'T 



•' z 



dz 



which can be rewritten as 



(3.32) 



B = 1 2 h b'w'(-h) - 1/2 h u’''b* (3.33) 

where if"b’''' is the downward surface flux of buoyancy : 

bnr=' = -“^(O) + ag pgCp RQJl + e' ^ ^ . 2'Yhj . (3.34) 

The radiation absorption Q(z) has been modeled as 

Q = yRQ/f^ (3.35) 

where y is the extinction coefficient for net solar radiation, and RQ^ is the short wave 
fraction of the net solar irradiance. This model assumes that the net long-wave solar 
radiation is absorbed at the surface, which leads to : 



-b'w'(O) = ag Pj^Cp [(1-R)Q3 - (Q,+ Qb + Qh)l + PgS(P-E) (3.36) 
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where (l-R)Qj corresponds to the short-wave incoming radiation. On the other hand, 
the short-wave radiation penetrates below the surface where it is absorbed following 
the exponential decay (3.35). 

c. Viscous dissipation 

For a fully turbulent mixed layer, viscous dissipation of the turbulence 
occurs primarily in the small eddies which are locally isotropic. The net rate of 
dissipation is modelled as follows; 



rO - , 

D = 1, „ c dz = m, < E> ' 
•-h-o ^ 



+ m^fh < E > 






or 



D = m, < E> ^' -(1 + Rq*' mj m, ir'/ < E> (3.3S) 

where Rg = u'-’/hf is a Rossby number for the turbulent boundary layer. The first 

term in D arises from the fact that the time scale of the largest eddies can be 

— 1 

proportional to the mixed layer depth divided by the rms turbulent velocity <E> 
The second one comes from the assumption that, in deeper boundary layers, planetarv’ 
rotation (time scale l;f) turns the mean shear direction with depth and thus influences 
the geometrical aspects of the integral scale. 
d. Redistribution of TKE 

The vertical integral of the pressure redistribution term 



,-0 p'(3u. 



= I .s V 



h-6 pg^X; 



az 



(3.39) 



is an important source or sink term for the individual TKE budgets, even though 
Rj R.> 4 - Rj = 0. The bulk formulation is 



R; = m,<E> (<E> - 3<U;'->) . 



(3.40) 



The rotational redistribution terms are assumed to be of higher order, and are 
neglected in this study. 
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5. Closure hypothesis 

Garwood [Ref. 1] achieves closure of the problem by formulating the following 
entrainment equation : 

ch dt = m_j<w'“> <E> ' hAB . (3.41) 



C. SUMMARY OF MODELED EQUATIONS 

We have six variables h, < B> , < U> , < V> , < w'-> , < q-> and a final set of 
six equations : 

• mean momentum equations ; 



h5<U>,5t + AAU 5h;^t = -fh<V> - u'w'(O) 



(3.42) 



h5<V> 5t + AAV^h/^t = -fh<U> - v'w'(O) 



(3.43) 



• mean buoyancy equation : 



h 5<B> ct + AAB 5h;5t = agQ^ ; PQCp - b'w'(O) 



(3.44) 



* horizontal integrated TKE component : 



5 , .3 (AU ‘ T- (AV ‘ dh 

--x-(h<q->) = m,u'’^ + ^ ^ 

2 ot ^ 3 2 at 



(3.45) 



- m2( < E> - 3 < w'"> ) < E> - 2m^ '3 ( < E> + m^'mj fh) < E> 
• vertical integrated TKE component : 



— 4-(h<w'^>) = 1/2 h b'w'(-h) - 1/2 h u*b* + m.,(<E> - 3 <vs''^> )< E> ‘ (3.46) 

2 Cl 



- m^ 3 (<E> + m.'mj fh)<E: 
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• entrainment equation : 



ch ct = < \v'" > 



-<E> ! hAB 



(3.47) 



The difi'erent constants used in the model are the following : 

• Representative density of air 

• Pq Representative density of sea-water 

• a.p CoefTicients of thermal expansion 

• g Gravitational acceleration 

• Cp Specific heat at constant pressure 

• Drag coefficient 

• f Coriolis parameter depending on latitude 

We can tune the model by specifying the dilTerent parameters m^, m_j, m^. m^ as well 
as the extinction coefficient for net solar radiation y, and the short-wave fraction of the 
solar radiation R. 

Initial profiles of temperature T(z) (°C vs cm) and salinity S{z) (g, kg vs cm) have 
to be provided. The depth increment in the model is 100 cm. Finally, as the time step 
is one hour, the following boundary conditions have to be prescribed hourly : 

• t Wind stress in dynes cm" 

• Qj Incoming solar radiation in W m" 

• Back radiation in W, m" 

• O Latent heat flux 

• Sensible heat fiux 

• E Evaporation in mm/hour 

• P Precipitation in mm/ hour 

The values of u* and u*b* can be computed using equations (3.17), (3.18), (3.30), 
(3.37) and (3.36). 
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IV. COUPLING THE TWO MODELS 



A. GENERALITIES 

The goal of this research is to couple the Oceanic Boundary Layer (OBL) model 
derived by Garwood in 1977 with an acoustic ray tracing program allowing us to 
analyze some effects of the atmospheric factors on the oceanic environment and, 
therefore, acoustic propagation in the ocean. 

The two previous chapters gave us a theoretical approach concerning these two 
models. Now, we are going to highlight the coupling of these two models, the inputs 
that the coupled model needs, and the outputs we can get from it. The block diagram 
of Figure 4.1 depicts the over plan, as followed in the ne.xt sections. 

Initially, we enter the OBL model with a set of boundary conditions, with some 
initial conditions, and with a set of parameters including the time-step of integration of 
the OBL model and the time interval between each resulting ray trace. From the 
predicted output of temperature and salinity profiles, we compute a sound-velocity 
profile which allows us to trace the paths of acoustic rays according to a set of 
geometrical parameters such as the source depth, the ma.ximum range, and the initial 
angle of a ray.' 

B. INITIAL CONDITIONS 

We have to initialize the coupled model with a temperature profile, that is, 
temperature T in Celsius versus depth z with a one-meter increment and a salinity 
profile, that is, salinity S in g kg versus depth z with a one-meter increment. These 
profiles can be obtained from climatological data [Ref 7], as the ones we used to 
analyze some applications of this coupled model, or from an XBT (e.xpendable 
bathythermograph) and a salinity profile from climatology. W'e can also assume a 
constant salinity profile based on an average salinity in the area studied because of the 
small effect of salinity variations on sound- velocity computation. 

"In all of the previous chapters, we wanted to be consistent with the coordinate 
systems used in each reference, that is, [Ref 4] for Chapter II and [Refs. 1.5] for 
Chapter III. In Chapter II, y is the downward vertical axis and, in Chapter 111, z is 
the upward vertical axis. Since the coupling is done through the output T and S 
profiles of the OBL model feeding the input profiles of the acoustic model, this 
notation conflict does not affect the coupling by itself nor the codes of the coupled 
programs. 
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4A 



C BOUNDARY CONDITIONS 



'I'he boundary’ conditions or atmospheric forcing factors we have to provide are 
listed below. 

The hourly solar radiation flux has to be given in W/nr versus time in hours. 
The climatology gives us generally a daily average In our study we will simulate 

SuVg 

a diurnal cycle for based on the following formula : 

Qj = Q^^ax if 0<t<12hours, (4.1) 

Qj = 0 if 12 < t < 24 hours. 

as depicted in Figure 4.2. 




Figure 4.2 Simulated diurnal cycle. 

A straight-forward integration leads to the following relation : 

Q = 7t Q . (4.2) 

^smax ^savg ^ ' 

Simulation of clouds is done by allowing only a part of this as input solar radiation 
(20 to 50% of Qj for example) at the sea surface. 
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The values of the latent flux of heat Q^, the sensible flux of heat Qj^ and the back 
radiation have to be provided on an hourly basis. In the simulations of the next 

chapter, we will pick some average values of the sum Qe'^Qh'''Qb 
climatology, depending on what time of the year we are working with. 

The evaporation rate EV in mm hour is also evaluated from climatological data, 
and precipitation rate PR in mm hour can be simulated. Hea\7 rain cases will be 
considered also because of some interesting effects. 



TABLE 2 

WIND SPEED AND WIND STRESS CORRESPONDENCE 



Wind Speed 




Wind Stress 


knots 


m's 


g/cm'sec^ 


5 


2.5 


0.106 


10 


5 


0.423 


15 


7.5 


0.959 


20 


10 


1.694 


25 


12.5 


2.647 


30 


15 


3.811 


35 


17.5 


5.188 


40 


20 


6.776 



The wind stress t in g cm sec^ has to be provided hourly also. In our 
simulations, we will consider both constant strong wind as well as light wind 
conditions. Based on equation (4.3) 

t = Pa Ccj L’jq" (4.3) 

3 

where an air density pa of 1.21 kg/m’’ and a drag coefficient C^ of 1.4x10' have been 
chosen. The correspondence between wind speed and wind stress is given in Table 2. 
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D. OTHER INPUT PARAMETERS 



The other parameters have to input in the coupled model arc the latitude of 
the area studied, the number of days NDA'VS during which we want to integrate the 
model, the interval of integration used by the model (1DIFF= Ihoiir in this study), and 
the frequency TI in hours with which we want to output some graphic results from the 
acoustic ray tracing program. 

E. OUTPUTS OF THE OBL MODEL 

From the OBL model we can get the hourly values of the predicted depth, 
temperature and salinity of the mi.xcd layer, and wc are able to follow the evolution of 
these dificrent parameters on tabulated outputs. 

The link with the acoustic model is done through the hourly computed profiles 
T{z) and S(z). Depending on the frequency TI, these profiles (renamed T(y) and S(y) 
for notational consistency as mentioned in a previous note) arc used to compute the 
sound-speed profile c(y) by a special subroutine SVEL. 

F. SOUND-SPEED COMPUTATION 

In the subroutine computing the sound-speed, we used the equation derived by 
Chen and Millero [Ref. S]. The sound-speed obtained agrees with the data of Del 
Grosso at one atmosphere [Ref. 9] and with the data of Wilson at high pressures 
[Ref. lOj. Other formulas for sound-speed are mentioned by Urick [Ref 11]. but the 
one wc used is most often suggested for oceanographic calculations. Further more, 
comparisons of graphic outputs do not reveal any significant dificrences between the 
results obtained by using these various equations. 

G. RAY TRACING 

The computed sound-speed profile c(y) is used as a main input in the acoustic 
model. Depending on the case we want to study, a set of geometrical parameters has 
to be provided such as the source depth YO, the maximum range lUANGE, and the 
maximum depth N. The bundle of rays to be plotted is defined by the number of rays 
M, the initial angle BO of the upper ray of the bundle, and the angle increment DELB 
between wo rays. Different equations for defining a bundle of rays can be set inside 
the subroutine. Thus, it is always possible to pick up a specific bundle of rays to 
highlight some effects of the atmospheric forcing on some precise type of rays and 
demonstrate some acoustic energy redistribution. 
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Since, in ihis research, we want to give an example of the application of the use 
of this coupled model, we will choose generally a ray-bundle width of 5 to 8 degrees 
which should correspond to an active sonar operating at a frequency in the range of 5 



Finally, the main output of the coupled model is a plot of the sound-speed profile 
aside an acoustic ray tracing for a given set of parameters defined previously. 

H. NOTE ON THE ALGORITHM USED IN THE ACOUSTIC MODEL 

Usually, in most acoustic models, the input sound-speed profile is approximated 
by a piecewise linear curve as depicted in Figure 4.3. 



kHz. 




C(y) 




Figure 4.3 Piecewise linear approximation 
of a typical sound-speed profile. 



Thus, the profile can be modeled mathematically by the following equation : 



c(y) = + gj(y - yj.j) 




(4.4) 
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and is composed of several distinct straight line segments, each with its own constant 
gradient gj. In each such delimited slab i, because of the constant gradient gj. the ray 
path is an arc of a circle for which an equation is easy to compute and fast to plot. 
Since a profile is usually divided into a few discrete slabs (for example, 10), a more 
computationally eiricient subroutine could be derived from the one discussed in 
Chapter II. This subroutine would invoke gradient sign checking and initial angle 
comparison to the horizontal only when a ray crosses a boundary between two slabs. 

In our study, this checking and comparison are done at each depth increment (1 
m) as well as the plotting of a new point belonging to the ray path. It is certainly a 
much slower algorithm than the previous one, but we had to consider that, even if the 
output profiles T(y) and S(y) of the Garwood OBL model can be approximated as 
piecewise linear, the output sound-speed profile c(y) cannot be considered as piecewise 
linear. This is because c(y) is a highly non-linear function of T. S, and y. This 
characteristic can be demonstrated by using a simple equation like the one derived by 
Coppens [Ref 1 2] 

c(T,S,y) = 1449.05 + 4.57T - 0.0521T“ + 0.00023T^ (4.5) 

+ (1.333 - 0.0126T)(S - 35) + 0.0163y . 

This may be rewritten as 

c = a + bT ^ cT- + dT^ + eS + ffS + gy . (4.6) 

Differentiation of c versus depth y yields 

dc dy = g + (b + 2cT + 3dT^+fS) dT/dy + (e+fS) dS dy . (4.7) 

Hence, in a particular slab i, even if (dT, dy)j and (dS'dy)j are constant, it can be seen 
from (4.7) that (dc;dy)j is not a constant because T and S are not constant in the slab, 
but are linear functions of depth y. 

Finally, the acoustic model developed in Chapter 11 permits the use of XBT 
profiles discretized with a depth increment of one meter or more. Use of climatology is 
often unrealistic due to multiple averaging and too coarse vertical resolution. 
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V. APPLICATION OF THE COUPLED MODEL TO A SPECIFIC AREA 



A. GENERALITIES 

In order to investigate the typical results we can obtain from this coupled-model, 
we shall simulate some cases in the Western Mediterranean Sea. The chart presented 
in Figure 5.1 shows the site to be studied: (p = 42.5‘’N, X = 007.5'’E. 




Figure 5.1 The northern region of the Algero-Provencal Basin. 

The Mediterranean Sea is one of the major, deep closed seas of the world. The 
water mass interacts with the open ocean only through the Strait of Gibraltar, which 
has a sill depth (300 m) much shallower than the average depth of the rest of these 
closed waters. The water below this sill remains almost homogeneous in temperature 
and salinity. The western Mediterranean is an area having very specific oceanographic 
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and atmospheric features. Figure 5.1 depicts the north part of the .*\lgero-Provencal 
Basin, which is one of the two large deep basins of the Western .Mediterranean. The 
basin bottom contains wide, flat abyssal plains, averaging 2700 meters in depth. The 
water mass is clearly bounded by land and islands, and this sea is well ventilated due to 
the quantity of deep water formed at the surface by winter cooling. As described in 
[Ref 13], during winter, the Alpine relief brings intrusions of polar continental air 
masses into contact with the surface waters of the northern part of the Mediterranean. 
Strong, cold, dry, continental winds (mistral) can blow for a few days, cooling the 
surface water and leading to a vigorous convective mixing. Such a strong sinking 
motion of surface water can mix and cool uniformly as much as the upper thousand 
meters of the water column. The sea-state is often characterued by a short wavelength 
swell which causes a very’ rough sea during storms. Another interesting result of the 
interaction between the atmosphere and the sea is that the Mediterranean water is 
characterized by high temperatures (13°C at 1000m) and high salinities (38.4-38.45 
g kg) which are surpassed only in the Red Sea. This is due to the large heat input and 
an excess of evaporation over fresh water input (precipitation and river contribution). 




Figure 5.2 Seasonal variation of the temperature profile 
in the northern part of the Western .Mediterranean. 

Figure 5.2 from (Ref 14] gives an example of the seasonal variations observed in the 
Ligurian Sea, that is, the northern part of the Algero-Provencal Basin. It can be 
observed that, in this region, the water column is completely isothermal in winter. A 
warmer surface layer evolves progressively during springtime. In autumn, it decreases 
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in temperature but increases in depth until winter conditions are reached again. The 
summer temperature profile is characterized by a top layer about 20 to 40 meters thick 
overlying a thermocline with large temperature gradients (up to 1 to 2^C, m ). 

The Mediterranean Sea is also of interest for the numerous naval operations 
conducted in that area. For this reason, it is of interest to study some of the features 
of the acoustic propagation that are related to atmospheric factors. For example, in 
this region, the heating by the sun of the upper layers of water, together with an 
absence of mixing by the wind, causes a strong near-surface negative sound-speed 
gradient to develop during the spring and summer months. This thermocline overlies 
isothermal water at greater depths. The result is a strong shallow internal sound 
channel (SOF.A.R channel: sound fixing and ranging) with its axial depth near 150 
meters. During the autumn, the profile returns to its winter time conditions, with 
isothermal water and positive sound-speed gradient extending to the surface of the 
water column, leading to a half-duct type of propagation. During the summer season 
the near-surface negative gradient and the resulting strong downward refraction greatly 
limit the detection range of surface ship hull-mounted sonars. By way of 
compensation, the summer time channel in the VIediterranean produces convergence 
zones for a near-surface source similar to the deep-ocean sound channel, although the 
range of the zone is much less because of the smaller vertical extent of the channel 
compared to the open ocean. Therefore, from the point of view of sound propagation, 
the important factor here is the existence of a predominantly isothermal and isohaline 
water mass at depth, which results in a sound-speed profile characterized by a deep 
constant-gradient section. Only the upper section of the profile will contain 
irregularities as depicted in Figures 5.3 and 5.4. See [Refs. 11,14] for more details. 

B. SUMMARY OF THE DIFFERENT PARAMETERS USED IN THIS STUDY 

In the simulations discussed in the next sections, we used the parameters listed in 
Table 3. Four different cases are to be analyzed, corresponding to the months of 
December, Februarv', June, and September. The boundar}’ conditions are assumed to 
be climatological and were obtained from [Refs. 16,17,18,19]. .A precipitation rate of 2 
mm hour over 12 hours will be assumed for hea\y rain cases. Values of 20 to 50% of 
0 will be used to simulate the reduction of radiation incident on the sea surface due to 
the presence of clouds. The initial profiles of temperature and salinity come from the 
climatology used by the Fleet Numerical Oceanographic Center (FNOC). Finally, we 
will avoid the complication of bottom reflection by assuming an infinitely deep ocean. 



SOUND VELOCITY (M/S) 



1500.0 1510.0 1520.0 1530.0 1540.0 




Figure 5.3 SSP for the Ligurian Sea, Western Mediterranean 
February to July from left to right. 
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SOUND VELOCITY (M/S) 



1500.0 1510.0 1520.0 1530.0 1540.0 




Figure 5.4 SSP for the Ligurian Sea, Western Mediterranean 
August to January from right to left. 
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C. DECEMBER CASE 



1. Source at 10 meters 

This case simulates a surface-ship hull-mounted sonar and we will consider 
strong wind (30 knots) as well as light wind (5 knots) conditions. 



TABLE 3 

MEAT FLUXES CLIMATOLOGY FOR THE .MED1TERR.\NEAN 

(W .M^) 





December 


February 


June 


September 


Qs 


130 


190 


250 


190 


Q 

^smax 


400 


600 


7S0 


600 


Qh-Qe 


190 


120 


80 


120 


Qb + Qh + Qe 


2S0 


190 


130 


190 


^nel 


-150 


0 


120 


0 


E\'(mm,'hr) 


0.3 


0.2 


0.1 


0.2 



a. Strong wind conditions {SO knots) 

The evolution of the mixed layer depth h for different cases is shown in 
Table 4. In the case of no clouds and no rain, we find that the mixed layer depth 
increases significantly due to the combined effects of the strong winter time surface 
cooling -150W m*) plus the effect of the wind stirring. 

Figure 5.5 shows the effect of this mixing and cooling on the acoustic ray 
paths. The first graph shows the acoustic propagation at time = 0 hour for a 
beamwidth 0 = 6®, and the second graph shows the propagation after 24 hours of 
strong wind. By counting the rays, we find that 39% of the acoustic energy has been 
redistributed. More rays are filling a deeper mi.xed layer, and the convergence zone 
{CZy has weakened. There is no effect on the CZ range which remains between IS 
and 20 km. 

Figure 5.6 shows, with more detail, the transformation of a bundle of 
refracted-surface reflected (R-SRi rays (initial angles between S7° and S9.3°) after 24 
hours of strong wind. For a total of 24 R-SR rays, 13 rays are now trapped in the 

^The term convergence zone (CZ) is applied to a phenomenon of focalization of 
the underwater sound in annular domains. 
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SSP (M/S) RANGE (KM) 

1507.0 15M.0 0.0 5.0 10.0 15.0 20.0 25.0 30.0 




Figure 5.5 December, strong wind, source at 10 meters: t = 0,24 hrs, 0 = 6^. 
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500.0 




Figure 5.6 December, transformation of R-SR rays, strong wind; t=0,24hrs. 
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mixed layer, increasing the possibility of detection on a shallow target for a surface- 
ship sonar. On the second graph of Figure 5.6, the weakening of the CZ is quite 
apparent. 



TABLE 4 

EVOLUTION OF THE MIXED LAYER DEPTH 
IN DECEMBER WITH A STRONG WIND 



no rain 
no clouds 

t(hr) h(m) 

0 37.0 

12 51.7 

24 60.3 

36 65.7 

48 70.8 

60 74.3 

72 78.2 



no ram 


no rain 


50% 


20% Qj 


h(m) 


h(m) 


37.0 


37.0 


52.5 


53.1 


61.2 


61.7 


67.0 


67.7 


72.0 


72.8 


76.0 


76.9 


79.5 


80.4 



heavy rain for 12 hrs 

20% Qj 

h(m) 

37.0 

50.6 
59.4 

65.7 
70.9 

75.1 

78.7 



From Table 4, we can sec that adding clouds to the simulation increases 
the surface cooling and thus accelerates a somew'hat the deepening of the mixed layer 
due to the added production of turbulence by wind stirring. The elTect of 
superimposing heavy rain for 12 hours is a slight damping of the mixed layer deepening 
rate. The main conclusion of these last test cases is that, with strong winds, the overall 
elTect of clouds and rain is not too significant and the general shape of the acoustic ray 
trace does not change. 

b. Light wind conditions (5 knots) 

Light wind conditions allow us to analyze the interesting elTect of heavy 
rain on the acoustic propagation. The evolution of the mixed layer depth is shown in 
Table 5 . 

Figures 5.7 and 5.8 show the elTect of a heavy rain on a winter time sound- 
speed profile. 

The gradient of the sound-speed versus depth can be expressed as : 
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Figure 5.7 December, light wind, heavy rain for 12hrs, 
source at 10 meters: t = 0,6 hrs, 6 = 5®. 



59 




Figure 5.8 December, light wind, heavy rain for I2hrs, 
source at 10 meters: t= 12,14 hrs, 0= 5®. 
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dc dy = dc 5T dT dy +oc dS dy + dc dp dp;dy 



(5.1) 



where the average values of the partial derivatives are the following : 

dc dT = 4.0 m, sec''C , (5.2) 

c*c,oS = 1.2 nrsec.g/kg , (5.3) 



and 

oc, i’p = 0.016 m sec/m . 



(5.4) 




The temperature term of the gradient dc/dy is dominant, and the salinity 
term is usually so small that it is very often neglected. In the case of hcavw rain, 
however, the term cc/oS dS, dy is not negligible and can be very important and even 
dominant. In the mixed layer, dT,dy = 0. Thus if heavw rain occurs, we get 
dS,dy> >0, and dc,'dy can be positive and much greater than usual. This yields a 
sharp "inversion" as depicted by the sound-speed profile of Figures 5.7 and 5.8. 
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Another way to examine the efTect of the rain is to consider the surface 
buoyancy flux Bq : 

^0 = «SQnet PCp - PgS(EV - PR) . (5.5) 

In winter time, the first term of is on average strongly negative. This corresponds 
to an upward buoyancy flux. During hea\w rain fall, EV-PR is negative, leading to a 
positive second term in Bq. This fact can be verified by investigating the last two 
columns of Table 4. Rain fall damps the upward buoyancy flux and, thus, decreases the 
mixed layer deepening. 

In Figures 5.7 and 5.S. we see the effect of heaw rain fall on the acoustic 
propagation at time 0. 6, 12, and 14 hours. In these graphs, the beamwidth has been 
set to 5°. After 6 hours of heavy rain, a high density of R-SR rays have been trapped 
in the mixed layer. After 12 hours, the effect is even much more important. Almost all 
the rays are trapped in a shallow mi.xed layer, and the convergence zone has 
disappeared. At that time, the only chance for a surface-ship sonar to get any 
detection is for a shallow target. 

However, this effect does not last veiy long after the rain stops. We can 
see in Table 5 that, in only two hours, the mixed layer depth almost returns to its 
initial value. We can also notice in Figures 5.7 and 5.8 that, only two hours after the 
end of the rain, the acoustic propagation is almost the same as for the initial case at 
t = 0. This is due to the strong winter time and night time surface cooling, even though 
the wind is light. This fact also demonstrates the strong effect of surface cooling on 
the mi.xed layer deepening in winter time. 

2. Source at 40 meters 

Now, we are going to study the case where the source is just below the initial 
mixed layer (37m). 

a. Strong n ind conditions {30 knots) 

In the case of no clouds and no rain. Figure 5.9 shows the evolution of the 
acoustic propagation over one day of strong wind for a beamwidth 0 of 6°. 

Initially, at time t = 0, all the rays are R-SR, and we have a perfect case of 
convergence zone (CZ) detection with a strong focusing effect at the smaller range of 
the CZ (around 15km). After 24 hours of strong wind, the deepening effect of the 
mixed layer redistributes 42% of the acoustic energy into a mixed layer type 
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Figure 5.9 December, strong wind, source at 40 meters: t=^0,24 hrs, 0 = 6°. 
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propagation. The CZ range stays the same, but the CZ is weakened and the focusing 
effect is greatly reduced. 

Figure 5.9 shows that the detection ability of a towed sonar at 40 m is 
much improved because of the deepening of the mi.xed layer mainly due to the wind 
stirring. For the purpose of highlighting this effect. Figure 5.10 extracts only the 
influenced ray bundle, that is. the bundle with a beamwidth of 2.5°. All the rays of this 
bundle are transformed from RR (refracted-refracted) type to mi.xed layer trapped. 

If we would have added clouds to the simulation, with or without rain, we 
would have found the same general shape for the acoustic ray trace. The fact confirms 
that, in winter time, the wind mixing has a dominating effect on the acoustic 
propagation. 

b. Light wind conditions (5 knots) 

For a source at 40 meters, even with light wind conditions, acoustic 
propagation is not very sensitive to hea\y rainfall. Figure 5.11 shows the ray trace at 
time 0 and after 12 hours of rain. The difference in shape between the two graphs is 
not ver\’ great. Thus, the effect of heavy' rain is to be considered significant, even with 
light wind conditions (which is not usually the case in winter time during storms 
passing by), for only a shallow transmitter, as seen in the previous section. 

3. Source at 100 meters 

If the source is sufficiently deep relative to the mixed layer depth, the 
atmospheric forcing has no effect on the acoustic propagation. Adding clouds and rain 
will not change the general shape of the acoustic ray paths. Figure 5.12 depicts the 
case of strong wind conditions with no clouds and no rain at t = 0 and 24 hours. For a 
beamwidth of 6°, the RR rays stay trapped in the shallow SOFAR channel and do not 
enter the mixed layer. This is because of the strong thermocline (strong negative 
gradient dc dy) underlying the mixed layer. For rays to enter the mixed layer, we 
would need steep initial angles at the source, which is not realistic for the usual sonar 
beamwidth. 

In Figure 5.13, with a source at 80 meters, an interesting effect on the ray 
paths is depicted. The rays almost appear to be reflected on the the sharp 
discontinuity created at the top of the thermocline because of the strong wind mixing. 

4. Source on the SOFAR axis (150 meters) 

Figure 5.14 shows the acoustic ray paths obtained when the depth of the 
source is set to 150 meters, that is. the optimum depth for the internal sound-channel 
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Figure 5.10 December, strong wind, source at 40 meters: t = 0,24 hrs, 0= 2.5°. 
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Figure 5.11 December, light wind, heavy rain, source at 40 meters; t = 0,12 hrs, 0 = 6°. 
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SSP (M/S) RRNGE (KM) 

1507.0 15M.0 0.0 5.0 10. 0 15.0 20.0 25.0 30.0 






(W) Hld3a 



(W) Hld3Q 



Figure 5.12 December, strong wind, source at 100 meters: t = 0,24 hrs, 0 = 6'’. 
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Figure 5.13 December, strong wind, source at 80 meters: t = 0,24 hrs, Pq = 91* to 93*. 
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Figure 5.14 December, strong wind, source at 150 m: t = 0,48 hrs, 0 = 9°. 
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propagation. Strong wind conditions were simulated for 48 hours and the beamwidth 
was set to 9®. The interesting acoustic propagation is due to the sharp positive and 
negative gradients dc/dy under and overlying the depth of the SOFAR axis. As 
mentioned previously, the source is deep enough for the acoustic propagation not to be 
influenced by the atmospheric forcing. A weak CZ appears at ranges between 15 to 18 
km, and these ranges are influenced almost not at all by the boundar}^ conditions. 

In that case, if we superimpose the second graphs of Figures 5.5 and 5.14, we 
can see the combination of the ray traces obtained with a source at 10 meters and a 
source at 150 meters. The water mass is almost completely filled with acoustic energy, 
leading to a high probability of detection of a target whose depth would be between 0 
and 300 meters. 

5. Source at 500 meters 

Figure 5.15 displays the case where the wind has been blowing for 24 hours 
with no clouds and no rain. We would have obtained the same general shape by 
varying the wind and (or) adding clouds and rain. As in the two previous cases, the 
source is deep enough for the acoustic propagation not to be influenced by the 
atmospheric factors. In any case, these graphs were provided to show the excellent 
shape of the acoustic ray trace in the case of a perfect surface reflection. 

6. Conclusions for December 

On short time scale, the atmospheric conditions have an efTect on the acoustic 
propagation only for a shallow source ( 10 to 40 meters in our simulations), that is, for 
a source in the initial mixed layer or just below. 

For such a shallow source, the wind has a dominating efTect over rain and 
clouds. Mainly, this effect is the transformation of R-SR rays to mixed layer trapped 
rays and the weakening of the convergence zone. 

The effect of the surface cooling is also important in winter time. The net 
upward heat flux was set to -150 W, nr in our simulations. This leads to an average 
deepening of the mixed layer untill the month of February when the entire water 
column is mixed. This was depicted in the light wind conditions cases, and it yields a 
half-duct type of acoustic propagation. 

The effect of heavy rain has been also depicted for light wind conditions and 
leads to a strong shallow inversion, trapping most of the rays emanating from a 
shallow source in a shallow surface duct. 
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Figure 5.15 December, strong wind, source at 500 m: t = 0,24 hrs, 0 = 6®. 
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D. FEBRUARY CASE 



1. Source at 10 meters 

a. Strong wind conditions {30 knots) 

Table 6 shows the evolution of the mixed layer depth for different cloud 
and rain conditions. In all the cases, the mixed layer deepens mainly because of the 
wind stirring and, as we saw previously in the December case, adding clouds to the 
simulation increases the deepening of the mixed layer, while the rain reduces slightly 
this cfiect. 



TABLE 6 

EVOLUTION OF THE MIXED LAYER DEPTH 
IX FEBRUARY WITH A STRONG WIN'D 





no rain 


no rain 


no rain 


hea\w rain for 12 hrs 




no clouds 


50% 


20% 


20% 


t(hr) 


h(m) 


h(m) 


h(m) 


h(m) 


0 


70.0 


70.0 


70.0 


70.0 


12 


77.6 


79.5 


80.5 


76.8 


24 


85.5 


87.3 


88.3 


84.9 


36 


89.1 


92.3 


94.1 


90.8 


48 


94.5 


97.7 


99.5 


96.3 


60 


96.9 


101.3 


103.8 


100.7 


72 


101.2 


105.6 


108.0 


104.7 



The initial mixed layer depth extracted by the OBL model is 70 meters. In 
this model, the criterion on which is based the choice of this depth is the following : 

|Ab| = IctgAT - PgAS| > 0.04 cm's“ (5.6) 

where the symbol A represents the variation of buoyancy, temperature, and salinity for 
an increment of one meter depth. A higher value for the criterion (0.05 for instance) 
would have increased all the values of the mixed layer depth, but the evolution would 
have had the same trend as the case depicted in Table 6. 
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However, independent of the choice of the criterion for the mixed layer 
depth, the sound-speed profile of Figure 5.16 shows that the water column is mixed 
almost from the top to the bottom with an average sound-speed gradient dc dy of 0.02 
sec'^ Thus, in February, the main characteristic of the acoustic propagation is a half- 
duct type of propagation. In the case illustrated by Figure 5.16, the beamwidth has 
been set to 12°. and the two acoustic ray traces show the situation at time 0 and 24 
hours. The wind does not affect the acoustic propagation since the water mass is 
already almost completly mixed at t = 0. Nearly identical ray path shapes would have 
been obtained by assuming clouds and rain. 

b. Light wind conditions (5-10 knots) 

Table 7. for light wind, no clouds, and no rain conditions, highlights the 
diurnal variability of the mi.xed layer depth and illustrates how the surface waters of the 
sea warm during the daytime hours on a sunny day, shallowing the mixed layer depth. 
At night, the deepening and cooling mixed layer returns to a state nearly the same as 
the initial one. These diurnal changes have a profound effect on sound transmission 
for a surface-ship sonar as depicted in Figures 5.17 and 5. IS. 

In Figure 5.17, we simulate a light wind of 5 knots with no clouds and no 
rain. The ray traces are shown at time instant 0 and 12 hours for a beamwidth of S°. 
At the end of the daytime (t= 12 hours), the general shape of the acoustic ray trace is 
almost the same as the one at the end of the night time (t = 0 hour), e.xcept for the 
ranges less than 6 km in shallow water. This is depicted in Figure 5.18 at t = S hours, 
that is, in the middle of the afternoon, for a beamwidth of 6°. This reduction of 
surface-ship echo ranging on a shallow target is a phenomenon often called the 
"afternoon elTect." 

Next, hea\T rain is included for 12 hours with light wind conditions. The 
effect is depicted at time 0 and 12 hours in Figure 5.19, where the wind speed was set 
equal 10 knots and the beamwidth to 6°. At the initial time, the acoustic energy 
propagates following a half-duct. At time 12 hours, we find that 62% of the rays are 
trapped in a shallow surface duct because of the strong discontinuity in the sound- 
speed profile due to the decrease of salinity in the upper layer of the sea caused by a 
weak wind stirring and damping of the turbulence by the strong downward surface 
buoyancy fiux associated with the rainfall. However, the overall shape of the acoustic 
ray trace is not changed significantly. There will be an improved possibility of 
detection for a shallow target. If the wind continues to blow at 10 knots, the effect 
remains even 12 hours after the rain has stopped. 
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Figure 5.16 February’, strong wind, source at 10 meters: t = 0,24 hrs, 6= 12°. 
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Figure 5.17 



Februar>', light wind, source at 10 meters: t = 0,12 hrs, 0 = 8®. 
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(W) HldGO 



(W) Hld3Q 



Figure 5.18 February, light wind, source at 10 meters: t = 0,8 hrs, 0 = 6°. 
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Figure 5.19 Februar>', light wind, hea \7 rain, source at 10 meters; t = 0,12 hrs, 0 = 8°. 
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2. Source at 200 meters 

In this case, we simulate the performance of a sonar towed at a depth of 200 
meters. In order to demonstrate that the atmospheric forcing has no effect on sound 
transmission from a deep source. Figure 5.20 shows the case with a heavy rain for 12 
hours together with a light wind of 10 knots. We can see that rain does not influence 
the acoustic propagation. Whatever the atmospheric conditions in February, a deep 
towed sonar is able to detect any target moving between 0 and 300 meters. Thus, the 
half-duct type of propagation is one of the best configurations for underwater detection 
by deep transmitters. 



TABLE 7 

EVOLUTION OF THE .MIXED LAYER DEPTH 
IN FEBRUARY WITH A LIGHT WIND 



wind 5 knots 




wind 5 knots 




wind 10 knots 




no rain 




heavy rain 




heavy rain 








for 12 hrs 




for 12 hrs 




no clouds 




20% 




20% 




t(hr) 


h(m) 


t(hr) 


h(m) 


t(hr) 


h(m) 


0 


70.0 


0 


70.0 


0 


70.0 


12 


7.6 


12 


2.1 


12 


12.9 


24 


69.0 


16 


3.3 


20 


14.0 


36 


7.6 


20 


7.9 


24 


16.8 


48 


69.1 


24 


67.8 


36 


33.7 


60 


7.6 


36 


70.9 


48 


70.8 


72 


69.3 


48 


73.2 


60 


73.0 



3. Conclusions for February 

As in December, the atmospheric conditions do not influence the acoustic 
propagation for rays emanating from a deep towed sonar. 

Compared to December, the SOFAR channel has disappeared, and there is a 
completely asymmetrical, surftce-bounded channel (or half-duct) inside which the 



78 




in 

\ 

IZ 



Q_ 

in 

in 






(W) Hld3Q 



(W) Hld3Q 



Figure 5.20 Februar}', light wind, hcavw rain, source at 200 meters: t = 0,12 hrs, 0=8^ 
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propagation is made of R-SR rays only. Because of this half-duct type of propagation, 
Februar}’ is one of the best months in the Mediterranean for underwater detection by 
botli deep and shallow transmitters. This generalisation not take in consideration the 
sea state which can generate noise and reflection loss when the sea is rough, because 
our acoustic model assumes a perfect surface reflection. 

With light wind conditions, surface-ship echo ranging can be poorer in the 
afternoon on a shallow target because of the "afternoon effect" due to surface water 
heating. 




Figure 5.21 Variation of sound-channel width d for similar temperature 
changes in (a) the Mediterranean and (b) the .^.tlantic. 



E. JUNE CASE 

Figures 5.3 and 5.4 show that the shape of the sound-speed profile stays 
approximatively the same during the entire suntmer season, that is, from June to 
September. Therefore, in this section, we will just simulate the atmospheric conditions 
and the acoustic propagation in June. All of our observations will be applicable to the 
other months of summer time even though the atmospheric factors change somewhat, 
as illustrated in Table 3. 

In summer, so much heat is added to the upper layer of the Mediterranean Sea 
that negative gradients in temperature and, thus, in sound-speed develop in the upper 



SO 



layer of the sea. This is underlain by isothermal water at a relatively shallow depth. 
The consequence is that a sound channel exists in summer in the .Mediterranean and is 
characterized as follows: 

The sound channel is, in all cases, strongly asymmetrical. The sound-speed 
profile contains a sharp discontinuity between sound-speed gradients of completely 
different orders of magnitude; 0.017 sec'^ on the lower side, due to the pressure effect, 
and -1 to -5 sec'* on the upper side in the thermocline. 



TABLE S 

EVOLUTION OF THE MIXED LAYER DEPTH 
IN JUNE WITH STRONG WINDS 





no rain 


no rain 


no rain 


heavy rain for 12 hrs 




no clouds 


50% 


20% 


207^0 Qj 


t(hr) 


h(m) 


h(m) 


h(m) 


h(m) 


0 


2.0 


2.0 


2.0 


2.0 


12 


20.3 


20.8 


21.1 


20.3 


24 


25.7 


26.2 


26.5 


25.8 


36 


29.4 


30.3 


30.8 


30.0 


48 


32.9 


33.9 


34.5 


33.6 


60 


35.5 


36.9 


37.7 


36.8 


72 


38.3 


39.8 


40.7 


39.8 



The sound channel axis is very' close to the surface, and the curvature of the 
sound-speed profile near this minimum is very great. 

The width of the sound channel is smaller than in the Atlantic and is extremely 
variable with the temperature of the surface water, as depicted in Figure 5.21. The 
usual width of the sound channel in summer is on the order of 1400 to 1800 meters. 
The depth of the Mediterranean is sufiicient for this channel not to be intercepted by 
the sea bottom in most cases. 
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1. Source at 10 meters 

a. Strong wind conditions (50 knots) 

The evolution of the mixed layer depth h for different cases is shown in 
Table 8. in the case of no clouds and no rain, we find that the mixed layer depth 
increases significantly during the first 12 hours because of the strong wind mixing in 
spite of the strong summer time heating. As seen in the previous sections, for short 
time scales, the wind is the dominating factor influencing both the mixed layer depth 
and the acoustic propagation. 

Figure 5.22 shows the acoustic ray trace at time 0 for a beamwidth 0 = 8°. 
After 24 hours of 30 knot winds, the mixed layer deepens to 25.7 meters. Even if the 
convergence zone weakens slightly, the range of the CZ is not affected. However, the 
acoustic propagation changes greatly for the less steep rays. Thirty percent of the rays 
are now trapped in a shallow mixed layer, improving slightly the probability of 
detection using a surface-ship sonar. However, the unfavorable conditions for 
underwater detection for a shallow source is well depicted in Figure 5.23. 



TABLE 9 

MAXIMUM AND MINIMUM .ML DEPTHS AND TEMPERATURES 
IN JUNE FOR LIGHT WIND CONDITIONS WITH NO CLOUDS 



day 2 


9a.m. 


h = 4.8m 


T=20.13°C 




5 p.m. 


h= 1.1m 


T=21.31°C 


day 3 


9a.m. 


h= 5.4m 


T=20.15°C 




5p.m. 


h= 1.1m 


T=21.33°C 


day 4 


9a.m. 


h= 5.6m 


T = 20.21°C 




5p.m. 


h = 1.1m 


T=21.39°C 



b. Light wind conditions (5 knots) 

Table 9 shows the evolution of the maximum and minimum mixed layer 
depths and high temperatures rise during light wind conditions with no clouds. The 
mixed layer responds to the diurnal cycle of solar radiation, shallowing in daytime 
because of solar heating, and deepening at night because of surface cooling. We notice 
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Figure 5.22 



June, strong wind, source at lOm: t = 0,24 hrs, 6 = 8®. 
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Figure 5.23 June, strong wind, source at 10m: t = 0,24 hrs, 0 = 8”. 
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that the general trend of the minimum and maximum temperature is to increase, 
leading progressively to the extreme sound-speed prolile of August {see I'igure 5.4). 
Xcthertheless, the mixed layer stays suniciently shallow not to influence the acoustic 
propagation, even for a shallow transmitter. 

•Adding clouds to the simulation (20% Q^) only deepens the mixed layer to 
9 meters after 48 hours, which has almost no elfect on the acoustic propagation. A 
longer period of' cloud covering would not be a realistic simulation for summertime in 
the Mediterranean. 

Finally, hea\w rain keeps the mixed layer shallow (around 1.5 meters), and 
sound propagation is not influenced. 

2. Source on the SOFAR a.xis (90 meters) 

Figure 5.24 simulates the case where the source is on the shallow sound- 
channel axis characterizing the Mediterranean summer time. The beamwidth is set to 
8°. The strong asymmetry of the sound channel is readily apparent. As previously 
mentioned, the wind does not influence the sound channel propagation for G = 8°. The 
rays at the source are not sufliciently steep to be affected and, even if the sound 
channel width decreases, the propagation does not change. 

However, in order to study the effect due to the decrease of the sound channel 
width. Figure 5.25 shows the ray traces for steep initial angles P = 8ff, 82'’, 83°. 97®. 
98°. and 99° at time 0 and after 48 hours of strong winds (30 knots). The channel 
width decreases from 1400 to 900 meters. The steepest rays (P=81°, 82°, 98°. and 99°) 
are transformed from RR type to R-SR type, but the influence of this phenomenon on 
the acoustic propagation is almost not perceptible with an acoustic ray tracing model 
as shown by Figure 5.25. A transmission loss model would be necessarx’ to account for 
the transmission loss due to the reflection on a rough sea surface. 

3. Conclusions for June 

In June, and for the reminder of the summer, the sound propagation for a 
deep transmitter is characterized by a sound-channel type of propagation and is not 
influenced by the atmospheric forcing on a short time scale of a few days. 

Because of the absence of a mixed layer and the presence of a sharp 
thermocline due to the important effect of surface heating, the detection capabilities of 
a surface-ship sonar are poor and limited to the convergence zone with ranges between 
30 and 35 km. However, strong winds can create a shallow mixed layer in which 
shallow targets could be detected as depicted in Figure 5.26, where rays emanating 
from both shallow and deep transmitters are traced. 
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Figure 5.24 June, strong wind, source at 90 meters: t = 0,24 hrs, 6 = 8*. 
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Figure 5.25 



June, strong wind, source at 90 meters, t = 0,48 hrs, P = 8r,82'’,83°,97‘’,98° 



87 






(W) Hld30 



(U) Hld30 



Figure 5.26 June, strong wind, sources at 10 and 90 meters 
t = 0,48 hrs, 0 = 8°. 
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VI. CONCLUSION 



On a short lime scale of a few hours to three days, the atmospheric conditions 
influence the acoustic propagation only for a shallow source. The ray trace emanating 
from a deep transmitter is almost not affected by the atmospheric forcing. This 
conclusion is independent of the type of propagation, depending on the time of the 
year (deep-sound channel, convergence zone, or half duct). 

The wind is the dominant atmospheric factor. For short time scales, the coupled 
model is not very sensitive to the other boundary conditions (solar irradiation, and 
rainfall). This is a favorable finding in the sense that the wind is the easiest parameter 
to measure and to forecast at sea. However, the importance of a hea\w rainfall has 
been demonstrated in the case of light wind conditions. 

The initial temperature profile plays a determinant role and has to be provided as 
accurately as possible. The coupled model is not very sensitive to the initial salinity 
profile. 

A. ADVANTAGES OF THE COUPLED MODEL 

The graphic output sequence from the coupled model shows qualitatively the 
acoustic propagation and its evolution under the forcing of specified atmospheric 
conditions, and an assumed initial temperature profile. 

The simplicity of the model makes it possible to be implemented on a simple 
desktop computer with some graphics capability. 

A current XBT, digitized by using a one-meter depth increment, can be used as 
input in the model. 

B. WEAKNESS OF THE COUPLED MODEL 

The fact that the coupled model does not account very well for the operating 
frequency is inherent to a ray tracing routine. The user of such a model will have to 
simulate a beamwidth at the source corresponding approximatively to a given 
frequency. 

The output of the coupled model does not provide such quantitative results as 
median detection range, convergence zone range, S N ratio at the receiver, and it does 
not account for scattering and attenuation. 
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Perfect surface reflection has been assumed whatever the sea-state. .An 
infinitively-deep ocean has been assumed in order to avoid the complication of bottom 
reflection. However, the model can handle perfect bottom reflection or absorption. 
Furthermore, the earth curvature has been neglected. 

The coupled model assumes horizontal homogeneity of the water mass, which we 
considered to be realistic for the horizontal ranges considered. 

Finally, we neglected horizontal advection, which is compatible with the short 
time scales studied, especially in the .Mediterranean where the currents arc weak. 

C. RECOMMENDATIONS 

Further simulations could be applied to different locations in the Atlantic. 
Pacific, or Indian Oceans, where climatology is more available. 

Real data could be used to initialize the coupled model and the boundary 
conditions could be computed from observed atmospheric and oceanic conditions. The 
resulting outputs of the model could be compared with corresponding simulations. 

Then, one might couple the OBL model with one of the different acoustic models 
used by the US Navy' (such as RAYVIODE. FACT, or PE models) in order to obtain 
quantitative results including MDR or CZR for usual operating frequencies and a 
variety of environmental conditions. 
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APPENDIX 

RAY TRACING SUBROUTINE USE EXAMPLE 



File giving an example of sound speed profile to be used with RAY : 



0. 1482. 

30. 1482.9 

55. 1484,1 

90. 1480.5 

150. 1482. 

225. 1484.1 

300. 1436,2 



C 



10 



C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



LCDR JACQUES FOURNIOL. FRENCH NAVY, 

PARAMETER (Y0=50 . ,B0=90 , ,DB= , 5 ,M=7 ,MM=6 ,N=300 ,NN=5000 , RAMGE=30 , ) 
DIMENSION BET(M) ,CC(0:MM) ,YY(0:MM) 

DIMENSION C(0:N) ,G(-1 :N) ,YC(0 :N) 

DIMENSION Y(0:NN) ,Z(0:NN) 

READ (05, 10) (YY(I) ,CC(I) ,I=0,MM) 

FORMAT(F7.2,2X,F7.2) 

CALL RAY(Y0,BET,B0,DB,M,CC,YY,MM,C,G,YC,N,Y,Z,NN, RANGE) 

STOP 

END 



SUBROUTINE USING THE MM+1 POINTS OF A STRAIGTH LINE SEGMENT 
CONTINUOUS PROFILE ( CC VS YY ) TO TRACE M ACOUSTIC RAYS ISSUED 
FROM A SOURCE AT A DEPTH YO , 

REF. "UNDERWATER ACOUSTICS-A LINEAR SYSTEMS THEORY APPROACH" 

BY ZIOMEK L. P237-238. ACADEMIC PRESS, ORLANDO, FLORIDA. 1985, 
PERFECT SURFACE AND BOTTOM REFLECTIONS. FLAT BOTTOM. 

BETAO : STARTING ANGLE OF RAY, REF VERTICAL, IN RAD. 

M : MB OF RAYS SENT, 

BET : ARRAY OF THE INITIAL ANGLES IN DEGREES AT THE SOURCE. 

BO : INITIAL ANGLE OF THE UPPER RAY OF THE BUNDLE. 

DB : INCREMENT OF INITIAL ANGLES IN BET. 

YO : SOURCE DEPTH IN METER. 

MM+1 : MB OF PTS IN PROVIDED SSP, 

CC : PROVIDED SSP IN M/SEC. 

YY : DEPTH OF THE PTS OF THIS SSP IN METERS. 

NN : INDEX OF RANGE, 

Y : DEPTH IN METERS AND 

Z : HOR. DISTANCE TRAVELLED BY A RAY IN KM, 

N ; MAX INTEGER VALUE OF DEPTH IN METER. 

C : SSP .1 METER INCREMENT . 
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C G ; GRADIENT OF SSP. 

C YC : ARRAY OF DEPTH FOR PLOTTING SSP. 

C RANGE : MAN. RANGE OF THE PLOTIN KM. 

C K : DEPTH INDEX. 

C N MUST BE THE INTEGER VALUE OF THE MAX. DEPTH YY(MM) OF THE 
C INITIAL SSP. 

C 

SUBROUTINE RAY(Y0 , BET , BO , DB ,M, CC , YY,MM, C ,G, YC ,N, Y,Z ,NN, RANGE) 

C 

DIMENSION BET(M) ,CC(0:MM) ,YY(0:MM) 

DIMENSION C(0 :N) ,G(-1 :N) ,YC(0 :N) 

DIMENSION Y(0:NN),Z(0:NN) 

DATA PI/3.141592654/ 

Y(0)=YY(0) 

YC(0)=YY(0) 

JMIN=0 

C(0)=CC(0) 

c 

C COMPUTE THE SSP C AND ITS GRADIENT G WITH 1 METER DEPTH INCREMENT 
DO 10 1=1, MM 

DYY=YY(I)-YY(I-1) 

JMAX=IFIX(YY(I)) 

GG=(CC(I)-CC(I-1))/DYY 
DO 30 J=JMIN, JMAX-1 
G(J)=GG 

C(J+1)=C(J)+GG 
30 YC(J+l)=FLOAT(J+l) 

JMIN=JMAX 
10 CONTINUE 
C 

C CREATE GRADIENT AT SFC AND BOTTOM TO HANDLE REFLECTION 
G(-1)=G(0) 

G(N)=G(N-1) 

C 

C DO 03 J=1,M 

C03 BET(J)=90.-FLOAT(J)*0.5 

DO 03 J=1,M 

03 BET(J)=BO-FLOAT(J)*DEL3 
C 

C RESEARCH OF CMIN AND CMAX FOR PLOTTING 
CMIN=1550. 

CMAX=1450. 

DO 40 1=0, N 
CMIM=AMIN1 ( CMIN , C ( I ) ) 

4 0 CMAX=AMAX 1 ( CMAX , C ( I ) ) 

CMIN=AINT(CMIN) 

CMAX=AINT(CMAX)+1. 
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c 

c 

c 

c 

c 

c 

c 



c 



c 

c 



c 

c 



PLOT THE SSP AND PREPARE PLOT FOR RAYS 
CALL COMPRS 
CALL CK41(4107) 

CALL TEK618 

CALL SHERPA( ' THESEOUT ' , ' B ' , 1 ) 

CALL HWROT( 'AUTO' ) 

CALL BLOWUP (0.745) 

CALL PAGE (11. ,5.5) 

CALL PHYSOR( 1.2,1.) 

CALL AREA2D(1. ,3.) 

CALL THKFRM(.Ol) 

CALL K.NAME('SSP (M/S)$ ' , -100) 

CALL YNAMS( 'DEPTH (M)$',100) 

CALL CROSS 
CALL YAXAMG(0.) 

CALL YTICKS(5) 

CALL XTICKS(IFIX(CMAX-CMIN)) 

CALL GRAF (CMIN,CNAX-CMIN,CMAX, FLOAT (N) ,-50. ,0. ) 
C.ALL CURVE (C,YC,N+1,0) 

CALL FRAME 
CALL ENDGR(O) 

CALL PHYSOR(2.7,l.) 

CALL AREA2D(7.3,3.) 

CALL XNAMEC RANGE (KM) $',-100) 

CALL YNAHE( '$' ,100) 

CALL YNONUM 
CALL CROSS 
C.ALL XTICKS(5) 

CALL GRAF(0. ,5. , RANGE, FLOAT (N) ,-50. ,0. ) 

DO 05 J=1,M 

Y(0)=Y0 

Z(0)=0. 

K=IFIX(Y0) 

BETA0=BET(J)/180.*PI 

C0=C(K) 

B=SIN(BETA0)/C0 



DO 20 1=1, NN 
NNN=I 

CHECK IF SFC REFLECTION 
IF (K.EQ.O) BETA0=PI-BETA0 



C 
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C CHECK IF BOTTOM REFLECTION 
IF (K.EQ.N) BETA0=PI-BETA0 
C 

IF ((G(K).GT.0).AND.(G(K-1).GT.0.)) THEN 
IF (BETAO.GT.PI/2.) THEN 
K=K-1 

BETA=PI-ASIN(B*C(K)) 

Y(I)=FLOAT(K) 

Z(I)=Z(I-1)+(COS(BETAO)-COS(BETA))/(B*G(K)) 

BETAO=BETA 

ELSE 

BETAC=ASIN(C(K)/C(K+1)) 

IF (BETAO.LT.BETAC) THEN 
K=K+1 

BETA=ASIN(B*C(K)) 

Y(I)=FLOAT(K) 

Z(I)=Z(I-1)+(COS(BETAO)-COS(BETA))/(B*G(K-1)) 

BETAO=BETA 

ELSE 

Y(I)=FLOAT(K) 

Z(I)=Z(I-1)+2>COS(BETAO)/(B*G(K)) 

BETA0=PI-BETA0 

ENDIF 

ENDIF 

ELSE 

IF {(G(K).LT.0.).AND.(G(K-1).LT.0.)) then 

IF (BETAO.LT.PI/2. ) THEN 
K=K+1 

BETA=ASIN(B*C(K)) 

Y(I)=FLOAT(K) 

Z( I)=Z(I-1)+ (COS (BETAO) -COS (BETA) )/(B*G(K-l)) 
BETAO=BETA 
ELSE 

BETAC=PI-ASIN(C(K)/C(K-1)) 

IF (BETAO.GT.BETAC) THEN 
K=K-1 

BETA=PI-ASIN(B*C(K) ) 

Y(I)=FLOAT(K) 

Z(I)=Z(I-1)+(COS(BETAO)-COS(BETA))/(B*G(K)) 

BETAO=BETA 

ELSE 

Y(I)=FLOAT(K) 

Z(I)=Z(I-1)+2.*COS(BETAO)/(B*G(K-1)) 

BETA0=PI-BETA0 

ENDIF 

ENDIF 

ELSE 
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IF((G(K) .LT.O. ) .AND. (G(K-l) .GT.O. ) ) THEN 
IF (BETAO.GT.PI/2.) THEN 
K=K-1 

BETA=PI-ASIN(B*C(K) ) 

Y(I)=FLOAT(K) 

Z(I)=Z(I-1)+(COS(BETAO)-COS(BETA))/(B*G(K) ) 
BETAO=BETA 
ELSE 
K=K+1 

BETA=ASIN(B*C(K)) 

Y(I)=FLOAT(K) 

Z(I)=Z(I-1)+(COS(BETAO)-COS(BETA))/(B*G(K-1)) 
BETAO=BETA 
END IF 
ELSE 

IF ((G(K). GT.O.). AND. (G(K-l). LT.O.)) THEN 
IF (EETAO.GT.PI/2. ) THEN 
BETAC=PI-ASIN(C(K)/C(K-1) ) 

IF (BETAO.GT.BETAC) THEN 
K=K-1 

BETA=PI-ASIN(B*C(K) ) 

Y(I)=FLOAT(K) 

Z(I )=Z (I-l )+(COS (BETAO) -COS (BETA) )/(B^G(K)) 
BETAO=BETA 
ELSE 

Y(I)=FLOAT(K) 

Z(I)=Z(I-1)+2.*COS(BET.AO)/(B*G(K-1)) 

BETAO=PI-BETAO 

ENDIF 

ELSE 

BETAC=ASIN(C(K)/C(K+1) ) 

IF (BETAO.lt. BETAC) THEN 
K=K+1 

BETA=ASIN(B*C(K)) 

Y(I)=FLOAT(K) 

Z(I)=Z(I-1)+(COS(BETAO) -COS (BETA) )/(B*G(K-l)) 
BETAO=BETA 
ELSE 

Y(I)=FLOAT(K) 

Z(I)=Z(I-1)+2.*COS(BETAO)/(B*G(K)) 

BETAO=PI-BETAO 

ENDIF 

ENDIF 

ENDIF 

ENDIF 

ENDIF 

ENDIF 
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20 



IF(Z(I) .GT.RANGE*1000. ) GO TO 21 
CONTINUE 
C 
C 

21 CONTINUE 

C 

C CONVERSION M TO KM 

DO 50 1=0, NNN 
50 Z(I)=Z(I)*0.001 

C 
C 

C PLOT THE RAYS 

CALL CURVE (Z,Y,NNN+1,0) 

C 

05 CONTINUE 

C 

CALL FRAME 

C CALL HEADIN( 'ACOUSTIC RAY TRACING$ ' , 100 , 2 . , 1 ) 

CALL ENDPL(O) 

CALL DONEPL 

STOP 

END 
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